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We present an ab initio based method which gives clear insight into the interplay between the hy- 
bridization, the coulomb exchange, and the crystal-field interactions, as the degree of 4/ localization 
is varied across a series of strongly correlated cerium systems. The results for the ordered magnetic 
moments, magnetic structure, and ordering temperatures are in excellent agreement with experi- 
ment, including the occurence of a moment collapse of non-Kondo origin. In contrast, standard 
ab initio density functional calculations fail to predict, even qualitatively, the trend of the unusual 
magentic properties. 
PACS: 71.27+a, 71.28+d, H.lOFd, 75.30Mb, 75.20Hr, 75.10Lp 



The difficulties and interest in treating strongly cor- 
related electron systems, and the consequences of cor- 
relation effects on magnetic behavior in the transitional 
4/ or 5/ localization regime, provide one of the central 
problems of condensed matter physics. 1-3 The transi- 
tional regime behavior is neither atomiclike nor itiner- 
ant. This gives rise to an extremely interesting range 
of phenomena, but also causes very great difficulties in 
treating the theory of these phenomena adequately, es- 
pecially in a way providing the ability to predict the be- 
havior of specific materials. 1-3 An adequate treatment 
requires treating the interelectronic coulomb interaction, 
i.e. the correlation effects, as constrained by exchange 
symmetry. 4-6 In this letter, we demonstrate an approach 
for treating these difficulties in predicting the interesting 
and complex behavior of an important series of cerium 
compounds. 

The isostructural (rock-salt structure) series of the 
cerium monopnictides CeX (X = P, As, Sb, Bi) and 
monochalcogenides (X=S, Se, Te) have become proto- 
type model systems for study, because of their unusual 
magnetic properties. 7-13 This series of strongly corre- 
lated electron systems offers the opportunity to vary sys- 
tematically, through chemical pressure, the lattice con- 
stant and the cerium-cerium separation on going down 
the pnictogen or chalcogen column, and hence tailor the 
degree of 4/ localization from the strongly correlated 
limit in the heavier systems to the weakly correlated 
limit in the lighter systems. 7-13 The calculated single- 
impurity Kondo temperature, Tk, presented below, is 
much smaller than the magnetic ordering temperature in 
these systems, and hence this series lies in the magnetic 
regime of the Kondo phase diagram. 14 Nevertheless, in 
this work we demonstrate that the sensitivity of the hy- 
bridization, coulomb exchange, and crystal-field interac- 
tions with the chemical environment gives rise to a vari- 
ety of unusual and interesting magnetic properties across 
the series, in agreement with experiment, including the 
occurence of a non-Kondo magnetic moment collapse. 



This class of cerium systems exhibits large magnetic 
anisotropy which changes from the < 001 > direction 
in the pnictides to the < 111 > direction in the chalco- 
genides. The low-temperature ordered magnetic moment 
increases with increasing lattice constant for the pnic- 
tides from 0.80/zb in CeP to 2.1fi B in CeSb and CeBi, 7-8 
while it decreases with increasing lattice constant for the 
chalcogenides from 0.57/j.b in CeS to 0.3/xb in CeTe. 7-9 
The magnetic moment collapse from CeSb to CeTe, with 
both systems having about the same lattice constant, 
is indicative of the sensitivity of the magnetic interac- 
tions to chemical environment. The experimentally ob- 
served low-temperature structure in CeBi and CeSb is 
the < 001 > antifcrromagnetic type IA (TTII), whereas 
in CeAs and CeP the structure is the < 001 > anti- 
ferromagnetic type I (TJ-)- 7 ' 15 The ordering temperature 
increases from 8K in CeP to 26K in CeBi for the pnic- 
tides, whereas it decreases from 8.4K in CeS to an un- 
usually low 2.2K in CeTe. 7-11 Another unusual feature 
of this series of cerium compounds is the large suppres- 
sion of the crystal field (CF) splitting of the Ce 3+ free- 
ion 4/5/2 multiplet from values expected from the be- 
havior of the heavier isostructural rare-earth pnictides 
or chalcogenides. 16 This can be understood 17 as arising 
from band-/ hybridization effects. In both the cerium 
monopnictides and monochalcogenides, the CF splitting 
between the T7 doublet and the Tg quartet decreases with 
increasing anion size, from 150 K for CeP to 10 K in CeBi 
and from 130 K for CeS to 30 K for CeTe, and it is about 
the same for the same row in both series, a rather sur- 
prising result in view of the additional valence electron 
on the chalcogen ion. 18 Neutron scattering experiments 
have shown 19 that the r7-doublet is the CF ground state 
in all the cerium pnictides and chalcogenides. 

In this paper we present material-predictive results 
from two ab initio based methods to study the change 
of magnetic properties across this series of cerium sys- 
tems. The first, ab initio based, method gives clear in- 
sight into the role of the three pertinent interactions: 1) 
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The band-/ hybridization- induced inter-cerium magnetic 
coupling; 2) the corresponding effects of band-/ coulomb 
exchange; and 3) the crystal-field interaction. This ap- 
proach allows us also to understand the interplay between 
these interactions as the degree of 4/ localization is varied 
across the series. The predictive calculations give results 
for the magnetic moments, magnetic structure, and or- 
dering temperatures in excellent agreement with experi- 
ment. Thus, this approach allows to understand and pre- 
dict a number of key features of observed behavior. First, 
is the very low moment and low ordering temperature of 
the antiferromagnetism observed in CeTe, an incipient 
heavy Fcrmion system. (For a review of theory and ex- 
perimental behavior of heavy Fermion systems see refer- 
ences 2,3,20,21.) This ab initio-based method, described 
below, predicts the magnetic moment and ordering tem- 
perature collapse from CeSb to CeTe, both systems hav- 
ing about the same lattice constant but CeTe having an 
additional p electron. The origin of the moment collapse 
is of non-Kondo origin. The earlier work of Shcng and 
Cooper 5 showed that this magnetic ordering reduction is 
accurately predicted without including any crystal-field 
effects. An erroneous statement appears in the recent re- 
view article by Santini et al. 22 stating that crystal-field 
effects played an important role in the calculated results 
of Sheng and Cooper. 5 This is incorrect, since crystal- 
field effects were not included in these calculations. We 
show in this paper that including the crystal-field effects 
modifies this behavior only quantitatively. Second, our 
results demonstrate that, while the band-/ coulomb ex- 
change mediated interatomic 4/4/ interactions dominate 
the magnetic behavior for the heavier systems, which 
are more localized because of the larger Ce-Ce sepa- 
ration, the opposite is true for the lighter, more delo- 
calized systems, where the hybridization-mediated cou- 
pling dominates the magnetic behavior. This reflects the 
great sensitivity of the relative importance of hybridiza- 
tion and coulomb exchange effects on magnetic ordering 
depending on the degree of 4/ localization. Third, we 
show that for the lighter more dclocalized systems the 
crystal-field interactions are much larger than the inter- 
ccrium interactions and hence dominate the magnetic be- 
havior. Finally, we predict the experimentally observed 
change of the ground-state magnetic structure from the 
< 001 > antiferromagnetic type IA (HID m CeBi and 
CeSb to the < 001 > antiferromagnetic type I (||) in 
CeAs and CeP. On the other hand, the second ab initio 
method, based on density functional theory within the 
local density approximation (LDA), 23 ' 24 fails to predict, 
even qualitatively, the trend of magnetic properties in 
this series of strongly correlated electron systems. 

The first, ab initio based, method employs the degener- 
ate Anderson lattice model which incorporates explicitly 
the hybridization and the coulomb exchange interactions 
on an equal footing 4,5 
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The parameters entering the model Hamiltonian, i.e., 
the band energies e^, the /-state energy e m , the on-site 
coulomb repulsion U, the hybridization matrix elements, 
Vkm, and the band-/ coulomb exchange J mm '(k, k') = 

(4>* k {nW m {r2) 7^ il) m >{n)<i>k'{ri)j are evaluated on a 
wholly ab initio basis from non-spin polarized full poten- 
tial linear muffin tin orbital 23 (FPLMTO) calculations. 
Here, r\i stands for \r\ — r 2 |; <j)k are the non-/basis states 
of the FPLMTO, and ip m are the localized / states. Be- 
cause the size of both the hybridization and coulomb ex- 
change matrix elements are much smaller (~ 0.1 eV) than 
the intraatomic coulomb interaction U (6cV), one can 
apply perturbation theory and evaluate the anisotropic 
two-ion 6X6 interaction matrices 4 ' 5 , £v™ 2 ™, 2 (R2 — Ri), 
which couple the two /-ions. The exchange interactions 
have three contributions: the wholly band-/ coulomb ex- 
change mediated interaction proportional to J^m' (k, k'), 
the wholly hybridization-mediated exchange interaction 
proportional to V£ m , and the cross term proportional to 
V k 2 m J mm '(k, k'). With the two- ion interactions having 
been determined, the low-temperature magnetic moment 
and the ordering temperature can be determined by use 
of a mean field calculation. 4 ' 5,8 We have previously ap- 
plied this ab initio based method to investigate the effect 
of hybridization-induced cerium-cerium interactions 4,17 
and the combined effect of both the hybridization and 
coulomb induced interactions 5 on the magnetic proper- 
ties of the heavier cerium pnictides and chalcogenides 
(CeBi, CeSb, and CeTe). However, these calculations 
did not take into account the crystal field interaction 
and employed a warped muffin-tin LMTO calculation for 
the parameters entering the model. The excellent agree- 
ment found 5 with experiment for the low-temperature 
magnetic moment and ordering temperature is relatively 
unaffected by the CF interaction, because the CF inter- 
action in the heavier cerium systems is smaller than the 
two- ion exchange interactions. 

The second method employs ab initio spin polarized 
electronic structure calculations based on the FPLMTO 
method 23 using 1) only spin polarization, with the or- 
bital polarization included only through the spin-orbit 
coupling, and 2) both the spin and orbital polariza- 
tion polarization. 24 In these calculations the 4/ states 
are treated as band states. The orbital polarization is 
taken into account by means of an eigenvalue shift 24 , 
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AV m = —E 3 L z mi, for the 4/ atom. Here, L z is the z- 
component of the cerium total orbital moment, m/ is the 
magnetic quantum number, and E 3 is the Racah param- 
eter evaluated self-consistently at each iteration. 

The crystalline field, which was neglected in the pre- 
vious calculations, 4 ' 5 is expected to affect the magnetic 
behavior considerably, if it is large. It is important to 
emphasize that since in the first method the 4/ states 
are treated as core states, they interact only with the 
spherical component of the effective one-electron poten- 
tial. Thus, the interaction of the atomic-like 4f state with 
the non-spherical components of the potential, giving rise 
to the CF splitting, Acf = <T 8 — £r 7 , is not included 
in the calculation of the model Hamiltonian parameters. 
In this paper, we generalize the first, ah initio based, 
method to include both the interatomic 4/-4/ coupling 
and the crystal-field interactions on an equal footing and 
to employ a full potential LMTO evaluation of the model 
Hamiltonian parameters. While the effect of the full po- 
tential on both the hybridization and coulomb exchange 
interactions is small, including the CF interaction will 
be shown to play a role as important as the interatomic 
4/-4/ interactions for understanding and predicting the 
overall trend in the unusual magnetic properties, as as 
one chemically tunes the degree of 4/ localization across 
this series of strongly correlated electron systems. The 
resultant Hamiltonian is 4 ' 5 

+B 4 J2(°1(i)+50t(i)) , 

(2) 

where the (0y ) are the two-ion 4/-4/ interaction ma- 
trices rotated to a common crystal-lattice axis, and the 
O4 and 0\ are the Stevens operators equivalents act- 
ing on the Ce 3+ free-ion 4/ 5 / 2 multiplet. 25 The CF split- 
ting is Acf — 360£?4; a positive B4 value gives the T7 
ground state, which is experimentally observed. 19 While 
our work in progress is aimed at evaluating the CF split- 
ting on a wholly ah initio basis, in the absence of an ah 
initio value of the CF interaction in this class of strongly 
correlated cerium systems, the A.cf is set to the experi- 
mental values listed in Table 3. 10 ' 19 

In Table |, we present the calculated values of the 
zero-temperature cerium magnetic moment from the 
FPLMTO electronic structure calculations. Listed in the 
table are values both with and without the orbital polar- 
ization correction taken into account. Note, the impor- 
tance of including the orbital polarization in these 4f cor- 
related electron systems. As expected, in all cases, the 
orbital polarization is found to be opposite to the spin 
polarization. Comparison of the total energies predicts 
that the magnetic anisotropy changes from the < 001 > 
direction in the pnictides to the < 111 > in the chalco- 
genides, in agreement with experiment. On the other 



TABLE I. Values of the calculated and experimental ~ 
low-temperature ordered magnetic moments for the cerium 
chalcogenides and pnictides in units of [s,b- Listed are the 
LMTO values for the spin moment (is, the orbital moment 
fiL, and total moment fj,, for the spin polarized only calcula- 
tion and for the calculation with spin polarization and orbital 
polarization correction. 
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-1.24 


1.99 


0.75 


0.57 
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-1.08 


1.02 


-0.06 


-1.26 


2.07 


0.81 


0.57 
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-1.15 


1.28 


0.07 


-1.31 


2.29 


0.98 


0.30 


CeP 


-0.80 


0.55 


-0.25 


-0.85 


1.27 


0.43 


0.80 


CeAs 


-0.84 


0.64 


-0.20 


-0.85 


1.42 


0.57 


0.80 


CeSb 


-0.86 


0.74 


-0.12 


-0.91 


1.61 


0.70 


2.06 


CeBi 


-0.86 


0.74 


-0.12 


-0.95 


1.69 


0.74 


2.10 



hand, except perhaps for the lighter chalcogenides (CeS 
and CeSe), comparison of the ab initio and experimen- 
tal values for the magnetic moment indicates the failure 
of the LDA calculations to treat properly the correlation 
effects of the 4/ states (treated as valence states) within 
the LDA as the degree of 4/ correlations increases in the 
heavier pnictide systems. Furthermore, the ab inito cal- 
culations fail to predict the large moment collapse from 
CeSb to CeTe, the latter being described as an incipient 
heavy Fermion system. 2 ' 3,20,21 

In Table ||, we list the values of the m = m' =1/2 ma- 
trix elements (characteristic matrix elements of the 6X6 
exchange interaction matrix) for the first three nearest- 
neighbor shells for the light (CeP and CeS) and the heav- 
ier compounds (CeSb and CeTe). Listed separately in 
this table are the three contributions to the interatomic 
4/-4/ interactions arising from band-/ hybridization (V 4 ), 
band-/ coulomb exchange (J 2 ), and the cross term. It is 
important to note that while the coulomb exchange me- 
diated interactions dominate the magnetic behavior for 
the heavier, more localized, 4/ systems, the opposite is 
true for the lighter, more delocalized, systems where the 
hybridization mediated interactions dominate the mag- 
netic behavior. This change of behavior of the inter- 
atomic 4/4/ interactions is a result of the sensitivity of 
the hybridization and coulomb exchange to the degree 
of 4/ localization. Equally important, is that while both 
first and second nearest-neighbor 4/4/ interactions are 
ferromagnetic for CeSb, there is an interplay between fer- 
romagnetic first nearest-neighbor and antiferromagnetic 
second nearest- neighbor interactions for CeTe. (These in- 
teractions are mediated via scattering of conduction elec- 
trons). This results in a saturated ordered moment for 
CeSb and in the ordered magnetic moment collapse for 
CeTe (see Table III). 

In order to determine whether the magnetic mo- 
ment collapse might be of Kondo origin, we have eval- 
uated the single-impurity Kondo temperature, 26 k^T^ 
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TABLE II. Values of the m = m! =1/2 matrix elements 
(characteristic matrix elements of the 6X6 interatomic 4/-4/ 

interaction matrix £; mi>m J (pt b — R a ), for the first, second, 

m a m a V ' > ' ' 

and third nearest-neighbor shells in degrees Kelvin. Listed 
are the values of the hybridization induced (E V i), cross terms 
(E V 2j), and pure coulomb exchange (-Ej 2 ) contributions. 
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-0.08 
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0.16 
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-0.16 


0.13 






CeSb CeTe 




E V 4 


E V 2J 
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E V 4 


E V 2J 


E.J2 


R= ( -2 5 o) 


0.70 


0.34 


7.30 
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-0.19 


2.90 


R = ( 1 ) 


2.07 


0.07 


10.21 


-0.19 


0.04 


-1.69 


R=(i||) 


-0.02 


-0.03 


0.40 


0.04 


-0.06 


-0.01 



TABLE III. Calculated values (from the first ab initio 
based method) of the zero-temperature ordered moment fj, in 
\xb, and the ordering temperature Tm in degrees Kelvin, with 
and without the crystalline field (CF) interaction across the 
cerium pnictide and chalcogenides series. Also listed are the 
experimental values 7-11 of fj, and Tat, and the CF splitting 19 
between the F7 ground state and the Fg state in degrees 
Kelvin. 
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CeS 
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0.73 


0.57 




1.0 


11.0 
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0.79 


0.57 
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14.0 


5.7 


CeTe 


32 


0.60 
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0.30 
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0.81 




14 
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0.85 
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37 
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1.80 


2.06 
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18 


17 


CeBi 
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2.10 


2.10 


2.10 




40 


40 


26 



= De 2 p< e f)| J( b f)i ) across the series. Here, D is the 
bandwidth of the conduction electron states, p(Ep) is 
the density of states of the conduction electrons at the 
Fermi energy, and J(Ep) is the conduction electron- 
/ exchange interaction at the Fermi energy, which has 
contributions both from the coulomb exchange interac- 
tion in Eq. (2), provided that it is negative, and the 
hybridization-induced exchange interaction \Jh v b{Ep)\ = 

{\E f -E^\)UEl-E F \+u) > where J hyb{E F ) < 0. We find 
that the coulomb exchange interaction in Eq. (2) eval- 
uated at the Fermi energy is positive across the entire 
series and hence cannot give rise to the Kondo effect. 
Thus, only the hybridization-induced exchange interac- 
tion, Jhyb(EF), can give rise to the Kondo effect. 26 Us- 
ing the ab initio values of the parameters entering the 
expression for T^, we find that Tk <C T rd across the 
entire series(Tx < 10 _4 K). The Kondo temperatures in 
the monopnictide series is smaller than that in the chaclo- 
genides, due to the fact that in the pnictides the Fermi en- 
ergy lies in the pseudogap, resulting in low p(Ep). These 
results suggest that the moment collapse from CeSb to 
CeTe is of non-Kondo origin. Rather, it results from an 
interplay of ferromagnetic and antiferromagnetic inter- 
atomic 4/-4/ interactions which arises purely from differ- 
ences in the underlying electronic structure. 

Listed in Table III are the calculated zero-temperature 
ordered moment and ordering temperature, Tn, from 
the first, ab initio based, method, with and without the 
CF interaction. It is clear that for the heavier systems 
(CeBi, CeSb, CeTe) the effect of the CF interaction on 
the magnetic moments is small, and it is slightly more 
pronounced on the ordering temperatures. This is due to 
the fact that for the more localized systems the CF in- 
teraction is smaller than the two-ion interactions. This is 
the reason that the previous calculations, 5 neglecting the 
CF interaction, gave results in very good agreement with 



experiment. On the other hand, for the lighter more delo- 
calized systems the CF interactions are much larger than 
the interatomic 4/-4/ interactions, and hence dominate 
the magnetic behavior. The overall decrease of the mag- 
netic moments in the presence of the CF interaction in all 
systems, arises from the mixing of the off-diagonal angu- 
lar momentum states | ±5/2 > and | ±3/2 > states from 
the CF interaction with ground state. Overall, we 
find that the first, ab initio based, approach which takes 
into account all three pertinent interactions (hybridiza- 
tion, coulomb exchange, and CF interactions) on an equal 
footing, yields results for both the zero-temperature mo- 
ment and the ordered temperature (a more stringent test 
for the theory) in excellent agreement with experiment. 

A final corroboration of the success of the first ab ini- 
tio based method is that it predicts the experimentally 
observed change of the ground-state magnetic structure 
from the < 001 > antiferromagnetic type IA (TTII) 
in CeBi and CeSb to the < 001 > antiferromagnetic 
type I (11) in CeAs and CeP. ' 15 More specifically, the 
sign [ferromagnetic (F) or antiferromagnetic (AF)] of the 
| ± 5/2 > matrix elements of the 6X6 exchange matrix 
determines the interplanar interaction between succes- 
sive (001) Ce planes. We find, that for the heavier com- 
pounds (CeBi and CeSb) the | ±5/2 > matrix elements 
of the coulomb exchange matrix are FM and hence favor 
the ttj], type, while in the lighter systems (CeAs and 
CeP) the | ± 5/2 > matrix elements of the hybridization- 
induced two-ion matrix are AF, and hence they favor the 

U type- 
In conclusion, we have applied two different, ab inito 
based and ab initio LDA, methods to study the dra- 
matic change of magnetic properties across a series of 
strongly correlated electron systems which offer the op- 
portunity to chemically tailor the different type of inter- 
actions (band-/ hybridization, band-/ coulomb exchange, 
and CF interactions), pertinent to the unusual magnetic 
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behavior. The first, ab initio based, approach which 
explicitly takes into account the interplay of the three 
pertinent interactions, gives results in excellent agree- 
ment with experiment for all compounds in the series, 
including the moment collapse from CeSb to CeTe and 
the trend of moments and ordering temperatures across 
the series. The remaining problem of determining on a 
wholly ab initio basis the suppressed crystal-field interac- 
tions in this class of systems poses a theoretical challenge 
for future theoretical work. On the other hand, the sec- 
ond, fully ab initio LDA, method gives good results for 
the lighter chalcogenide systems, but it entirely fails to 
give, even qualitatively, the trend of the unusual mag- 
netic behavior. 
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